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C. elegans is extensively used to study the Wnt-pathway and most of the core-signalling components are 
known. Four p-catenins are important gene expression regulators in Wnt- signalling. One of these, bar-1, is 
part of the canonical Wnt-pathway. Together with Wnt effector pop-1, bar-1 forms a transcription 
activation complex which regulates the transcription of downstream genes. The effects of bar-1 
loss-of- function mutations on many phenotypes have been studied well. However, the effects on global gene 
expression are unknown. Here we report the effects of a loss- of- function mutation bar-1 (gaSO). By analysing 
the transcriptome and developmental phenotyping we show that bar-l(ga80) impairs developmental 
timing. This developmental difference confounds the comparison of the gene expression profile between the 
mutant and the reference strain. When corrected for this difference it was possible to identify genes that 
were directly affected by the bar-1 mutation. We show that the Wnt-pathway itself is activated, as well as 
transcription factors e/t-3, pqm-1, mdl-1 and pha-4 and their associated genes. The outcomes imply that this 
response compensates for the loss of functional bar-1. Altogether we show that bar-1 loss-of function leads 
to delayed development possibly caused by an induction of a stress response, reflected by daf-16 activated 
genes. 

The Wnt/p-catenin pathway is highly conserved across metazoans and is essential for many cellular functions 
like cell specialization, cellular migration, adhesion and development. Although a key pathway in inverte- 
brates and vertebrates, the Wnt- signalling pathway yet remains to be fully elucidated^"^. A better under- 
standing is required not only from a fundamental biological point of view, but also because it can be important for 
developing drugs and medical treatments of Wnt associated diseases such as bone diseases and colorectal cancer 
(reviewed by^'^). 

The model worm Caenorhabditis elegans, a widely studied human model species, has the canonical Wnt- 
signalling pathway and a variation on this pathway, the asymmetrical cell division pathway. The Wnt pathway has 
five known Wnt genes: mom-2, cwn-1, cwn-2, lin-44 and egl-20^. In the canonical Wnt pathway, the cellular 
abundance of free P-catenin is controlled by a protein destruction complex which targets free p-catenin for 
proteasomal degradation. Activation of canonical Wnt- signalling, whereby Wnt binds to a Frizzled/LRP co- 
receptor, inactivates the destruction complex leading to accumulation of free P-catenin which then functions as a 
nuclear transcriptional activator. In C. elegans four distinct P-catenins have been identified. Bar-1 is part of the 
evolutionarily most conserved pathway^"^, whereas wrm-P°, hmp-2^\ and sys-P^ function in a variant of the Wnt 
pathway regulating asymmetrical cell divisions^'^. These P-catenins all seem to play distinct roles in the worm^l 

BAR-1 functions in the post-embryonic stage^'^^'^^ where it forms a transcription activation complex with the 
Wnt effector similar to the TCF/p-catenin complex in flies and vertebrates^^. BAR-1 is regulated by the 

axin-like protein PRY-1, the GSK3P homolog SGG-1 and the APC-like protein APR-1^^'^^. Among the processes 
influenced by BAR-1 are P12 cell fate specification ^^'^^ and mab-5 expression in the neuroblast QL^^'^°. Moreover, 
BAR-1 is involved in vulval precursor cell specification in the early LI stage of C. elegans through transcriptional 
activation of the Hox gene lin-39. In addition, BAR-1 plays a role in cell fate specification of the vulva during the 
early L3 staged The mutation used in this study, bar-l(ga80), affects vulval precursor cell induction which results 
in an incomplete vulva, a protruding vulva (pvl) and egg-laying defects (egl)^. 

Even though many developmental processes in which bar-1 is involved are known, the effect of bar-1 on gene 
expression is poorly understood. We studied gene expression patterns during the fourth larval stage (L4) of the 
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worm strain EW15 carrying the P-catenin-loss-of-function point 
mutation bar-l(ga80). This mutation causes a Glu to Stop codon 
change at amino acid 97 of the predicted BAR- 1 protein. By analysing 
the transcriptome and developmental phenotyping we show that 
bar-l(ga80) impairs developmental timing. Moreover we found that 
without a functioning bar-1 ~ 7,500 genes were affected. Our results 
further suggest that the loss of bar-1 is partially compensated by 
redundancy in the Wnt- signalling pathway, pointing towards a feed- 
back mechanism between p-catenin activation and the expression 
levels of Wnt- signalling pathway encoding genes. 

Methods 

Strains. The following strains were used: EW15 (Bristol N2 strain, carrying the 
mutation bar-l(ga80)) and wild type Bristol N2. Upon arrival in the lab EW15 was 
outcrossed at least 4 times. Worms were kept in maintenance at 12°C and before 
experiments started, populations were cleared of males and grown at 20°C until all 
worms were gravid. All experiments were conducted at 20°C. 

Microarray experiment. Strains were stage synchronized by bleaching^\ then grown 
on 9 cm NGM Petri dishes seeded with E. coli OP50. Worms were rinsed of the plates 
with M9 buffer 48 hours after synchronization, snap-frozen in liquid nitrogen and 
stored at — 80°C before performing gene expression profiling using microarray 
analysis. The experiment was performed in three independent duplicates. 

Egg laying experiment. Strains were synchronised by bleaching, after which 
approximately 200 eggs from N2 or 300 eggs from EW15 were transferred to a fresh 
9 cm NGM dish containing E. coli OP50. Starting at 58 hours after bleaching, the 
production of eggs by the adult worms was observed. Eggs on the plate were scored as: 
0 (no eggs), 1 (first 1-100 eggs, first worms have started laying), 2 (100-200 eggs, 
multiple worms are laying eggs), 3 (>200, most worms are laying eggs, first egg 
clusters appear) or 4 (many eggs and the first eggs are hatching) . The scoring was done 
every hour. The experiment was performed in three independent duplicates. 

Egg hatching experiment. Strains were synchronised by bleaching, after which as 
many eggs as possible (up to 700 eggs) were placed on a fresh 9 cm NGM dish 
containing E. coli OP50. Observations started immediately after bleaching, for every 
30 minutes. Once an egg hatched, the LI larvae was picked and counted. This was 
continued until all eggs hatched. The experiment was performed three times. 

Microarray sample preparation, scanning and normalization. mRNA isolation was 
performed using the RNeasy Micro Kit from Qiagen (Hilden, Germany), following 
the 'Purification of Total RNA from Animal and Human Tissues' -protocol provided 
with the kit. After this, the 'Two-Color Microarray-Based Gene Expression Analysis; 
Low Input Quick Amp Labeling' -protocol, version 6.0 from Agilent (Agilent 
Technologies, Santa Clara, CA, USA) was followed, starting from step 5. The 
microarrays used were C. elegans (V2) Gene Expression Microarray 4 X 44K slides, 
manufactured by Agilent. Input of total RNA was approximately 200 ng for each 
replicate. Three independent duplicates per strain were measured. The microarrays 
were scanned using an Agilent High Resolution C Scanner, using the settings as 
recommended in the above mentioned manual. Data was extracted with the Agilent 
Feature Extraction Software version 10.5, following manufacturers' guidelines. For 
normalization the Limma package for the "R" environment (version 2.13.1 x 64) was 
used. No background correction of the RNA-array data was performed as 
recommended by^^. For within-array normalization of the RNA-array data the Loess 
method was used and for between-array normalization the Quantile method was 
used. The obtained log2 normalized intensities (single channel data) were used for 
further analysis. 

Statistical analyses. All statistical analyses were performed using the statistical 
programming language "R" (version 2.13.1 x 64). A linear model was used to 
determine the effect and significance of the genotype on the expression levels (probe 
intensity ~ genotype + error). Using permutations of the original data in the same 
linear model, we determined thresholds adjusted for multiple testing (FDR 0.05: 
-loglO(p) > 2; FDR 0.01: -loglO(p) > 3). To correct for the developmental 
difference between bar-l(ga80) and N2 we used the developmental gene expression 
data from Snoek et al. (2014)^^ together with the gene expression data generated for 
this study {bar-l(ga80) vs. N2) in one linear model (probe intensity ~ sample age + 
genotype + error). The intensities were corrected for batch effect and for sample age 
we used an age of 44 hours for the bar-l(ga80) samples (as estimated), and 48 hours 
for the N2 samples. For the samples from Snoek et al. (2014) their original ages were 
used (44 to 58 hours after synchronisation). Genes with p > 0.05 for the "sample age" 
were selected as genes without a developmental effect. Genes with for which the 
"sample age" effect was opposite to the bar-1 effect were selected as genes with "effect 
which was opposite of what one would expect in a relatively slower developed bar- 
l(ga80) mutant". 

Enrichment tests were done using a hyper geometric test on the genes with a 
significant bar-l(ga80) effect (— loglO(p) > 2), excluding those with a developmental 
effect (unless stated otherwise). eQTL enrichments were done by selecting the genes 



with a significant (— loglO(p) > 3) linkage to each locus and comparing those against 
the genes affected by the mutant^"^. 

Datasets used. The GO -annotation, anatomy terms, protein domains and gene 
classes were obtained via Wormmart (www.caprica.caltech.edu:9002/biomart/ 
martview/) of the WS220 wormbase release. Genes from Wormbook chapters were 
obtained from the 2012 version of Wormbook (www.wormbook.org). Expression 
QTLs (eQTLs) were obtained from WormQTL (www.wormqtl.org;^^'^^ *^ using the 
data from^^"^''. Transcription factor binding sites were obtained from^" Binding sites 
from DAF-16 were obtained from modENCODE release #32 (www.modencode. 
org;^^), and mapped to transcription start sites according to Tepper et al. (2013)^^. 
KEGG pathways were obtained from Release 65.0 of the Kyoto Encyclopedia of Genes 
and Genomes (www.genome.jp/kegg/). 

Network visualization. The network of the transcription factors and their targets was 
visualized using Cytoscape (version 2.8.2)^^. 

Data storage. All data was stored in WormQTL (www.wormqtl.org)^^'^^'*^. 

Results and discussion 

bar-1 (gaSO) affects gene expression and slows development. We 

compared the transcriptomes of worms (age of 48 h) of N2 to the 
bar-1 (gaSO) mutant strain of the same age by microarray analysis. Of 
the 20,887 genes that were tested on the microarray, 5,772 genes were 
differentially expressed (-loglO(p) > 2.0 at FDR = 0.05). In bar- 
l(ga80), 51% was down-regulated (2,927 genes), and 49% was up- 
regulated (2,855 genes) compared to N2. During the initial analyses, 
we noticed that many of the differentially expressed genes were 
related to development, for example genes encoding for collagens 
and vitellogenins. Recently we reported that genome wide gene 
expression can rapidly and massively change during the L4 stage^^ 
To test for the developmental difference within the L4 stage, we 
compared the differentially expressed genes between N2 and bar- 
l(ga80) with the set of genes reported by Snoek et al. (2014)^^ 
(Figure 1). The differentially expressed genes between bar-l(ga80) 
and N2 were enriched for genes changing during L4 development 
(hypergeometric test, p < 1*10"^°°). To exactly pinpoint the 
developmental delay, we used the expression levels which have a 
linear correlation with L4 developmental timing to estimate the 
developmental age of the bar-l(ga80) and N2 samples. Even 
though all RNA samples from both genotypes were taken at 
48 hours after synchronisation we found that the bar-l(ga80) 
worms developed more slowly (44.4 h ± 0.9) than the N2 samples 
(47.7 h ± 0.8) (two-sided t-test, p = 6*10"') (Supplementary figure 
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Figure 1 | Genes affected by bar-1 and development. Volcano plot 
showing the effects and significance of the transcriptome comparison 
between bar-l(ga80) and N2. The black dots represent the spots on the 
array, the log2 effect between N2 and bar-1 (gaSO) is shown versus the LOD 
score. The yellow x indicate spots of genes affected by developmental 
effects during L4 development^^ 
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Figure 2 | Age-estimation and start of egg-laying in N2 and bar-l(ga80). 

A) Age estimates for N2 (blue) and bar-l(ga80) (yellow), based on 
linearly differentially expressed genes during L4 development^^. The bar- 
l(ga80) mutant is estimated significantly younger than N2 (p < 1*10"^). 

B) Start of egg-laying in hours after synchronization for N2 (blue) and 
bar-l(ga80). Again, the difference between the two strains is significant 
(p < 1*10"'). 

A and Figure 2 A) indicating a developmental delay of 3.3 hours after 
48 hours. 

To confirm this discrepancy we recorded the time until the first 
eggs were laid in N2 and bar-l(ga80). N2 started laying eggs at 
—62 hours and bar-l(ga80) started laying eggs at —68 hours (two- 
sided t-test, p = 4*10"^) (Figure 2B). To investigate whether the 
delay was caused in part by delayed hatching or slow embryonic 
development, the time until hatching after synchronization was 
determined. No difference was found between N2 and bar-l(ga80) 
in time from synchronising the eggs and hatching of those eggs 
(Supplementary figure B). 

The developmental delay of bar-l(ga80) increased over time (0 h at 
0 h, — 3.3 h at 48 h and 6 h at 62 hour). This implies that the muta- 
tion affected the entire developmental period from egg to adult. Our 
results show that bar-1 (gaSO) does not affect a single developmental 
stage, because than the developmental difference between bar-l(ga80) 
and N2 would remain constant during the subsequent stages. 

Analysis incorporating developmental effects. To exclude the 
effects of the developmental delay of bar-l(ga80) from other effects 
of bar-l(ga80) on gene expression, we included the transcriptional 
effects during L4 development in the analysis (Figure 3). Here we 
found 7,557 (FDR = 0.05) genes to be affected by the bar-l(ga80) 
mutation either with or without a development effect. Of these genes, 
3,920 were up-regulated and 3,637 were down -regulated in the bar- 
l(ga80) mutant (Supplementary figure C). 

As developmental effects were very strong and affected many 
genes^^'^'*, we selected those genes that did not have a developmental 
effect (P < 0.05) or an effect which was opposite of what one would 
expect in a relatively slower developed bar-l(ga80) mutant 
(Supplemental figure C). We also selected on effect size (>0.5 or 
< — 0.5) which resulted in 710 down- and 425 up-regulated genes 
compared to N2 (FDR = 0.05; Supplement Table 1). 

BAR-1 strongly affects collagens and hedge-hog signalling. The set 

of genes down- regulated in bar-l(ga80) compared to N2 (Supple- 
ment Table 1) contains many non-annotated genes. These genes 
could complement the genes with known functions, but could also 
constitute new functions. Furthermore genes like mai-1, dao-4,pho- 
11, sta-2, plc-2, pes-8, cnp-2, hmit-l.l, hmit-1.2, gcy-32, nlp-23 and 
fkb-5 have a strongly reduced expression in bar-l(ga80). These genes 
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Figure 3 | Effects on Wnt pathway genes. Members of the Wnt pathway 
are shown by their function. The first column shows the log2 -ratio between 
the bar-1 mutant and N2, blue indicates lower in bar-ly yellow higher in 
bar-L The second column shows the significance (-loglO(p)) of the bar-1 
effect. The third column shows the log2-ratio during L4 stage (44 to 
58 hours at 20°C)^^. Two red lines indicate the relative development of the 
bar-1 and N2 samples taken at 48 hours post hatching. Last column shows 
the significance (-loglO(p)) of the developmental effect. 

might function together with the strongly down-regulated groups of 
genes, like collagens {rol-1, bli-2, bli-1, dpy-3, lon-3), co/-type genes 
(col-175, -38, -71, -120, -40, -49, -138, -110, -97, -79, -70) or other 
cuticle related components {cutl-18, cutl-28, mlt-18, mltn-12, nas-27 
and gly-1). Components of hedge-hog (hh) signalling were also much 
lower expressed in bar-l(ga80) such as: the warthog genes wrt-6 and 
wrt-4; groundhog-like genes grl-15, grl-5 and grl-14; hedgehog-like 
genes, hog-1, grd-2, grd-1 and grd- 12. This shows that hh-signalling is 
affected by bar-l(ga80) mutation. Taken together, BAR-1 activity is 
most likely required for activation of collagens and other cuticle 
genes as well as genes involved in hedge-hog signalling. 

Expression of Wnt-signalling components. The expression of most 
Wnt-signalling components changed during development and were 
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affected by bar-l(ga80) mutation (Figure 3) (core-Wnt pathway 
genes selected by^). All four P-catenins were differentially 
expressed in bar-l(ga80). Expression of bar-1 was lower whereas 
wrm-1, hmp-1 and sys-1 all showed a slight increase in expression 
(see 4 upper blocks in the first column of Figure 3). In N2, these three 
P-catenins showed increased expression levels during L4 develop- 
ment whereas expression of bar-1 hardly changed throughout the L4 
stage in N2 (p = 0.042; see 4 upper blocks in the third column of 
Figure 3). 

Of the five Wnt genes in C. eleganSy only mom-2 was higher 
expressed in bar-l(ga80). For the other Wnt genes no significant 
effect was found. Of the Wnt genes, only mom-2 expression increased 
during L4 development, whereas expression of cwn-1, cwn-2 and lin- 
44 decreased (p < 0.001 in all cases). Egl-20 was not differentially 
expressed throughout L4 and was also not affected by the bar-1 
mutation. 

The members of the destruction complex showed no or only 
minor expression differences between bar-l(ga80) and N2 during 
development. One of the dishevelled genes, dsh-1, was affected by 
the bar-1 mutation. This gene was higher expressed in bar-l(ga80). It 
was also the only gene of the three dishevelled genes that did not 
show a change in expression. 

The Wnt receptors cfz-2, mig-1 and lin-18 {ryk/ derailed) were 
higher expressed in bar-l(ga80), and their expression did not change 
during development. The two other Wnt-receptors were slightly 
affected by development. Expression of mom-5 increased during 
L4 development whereas lin-17 expression decreased. Of the other 
genes, sfrp-l, an extracellular active Wnt- inhibit or ^^ had a lower 
expression in the bar-1 mutant compared to N2. The expression of 
sfrp-1 decreased during L4 development in N2. Even though bar- 
l(ga80) showed a developmental delay, the expression of sfrp-1 
decreased compared to N2. Some of the transcription factors known 
to be involved in the Wnt-pathway, like EOR-P^ and EGL-27'^ (or 
associated with the Wnt pathway, like SKN-P^), were up regulated, 
but their targets were not enriched for in the differentially expressed 
genes set (hypergeometric test, p > 0.1). 

To summarize, the up -regulation of the other three P-catenins in 
bar-1 (ga80) probably compensates for the loss of a functional BAR- 1 
(Supplement figure C). Intriguingly, not only the P-catenins were 
higher expressed, but also the Wnt-receptors. Furthermore, the sfrp- 
1 gene was down-regulated. Thus, a lack of bar-1 also affects the 
Wnt- signalling pathway upstream, which could point to a feed- 
back-mechanism. As the Wnt encoding gene mom-2, the Wnt-recep- 
tors cfz-2, mig-1 and lin-18 and the dishevelled gene dsh-1 were also 
up-regulated, our results imply that the Wnt- signalling pathway 
itself was activated following the knock-down of bar-1. 

Natural genetic variation in Wnt-pathway genes. All genes part of, 
or associated with, the Wnt-pathway (Figure 3) are polymorphic 
across many other C. elegans wild type strains^^"^\ Between the 
two most frequently studied wild types N2 and CB4856 these 
polymorphisms lead to an amino acid change in almost 50% of the 
proteins (Supplement text 1). Furthermore enrichments of 
expression Quantitative Trait Loci (eQTL) of genes with affected 
transcript levels by bar-l(ga80) suggest that polymorphic loci 
between CB4856 and N2 downstream of or modulated by bar-1 
and Wnt-signalling might be present (Supplement text 1). This 
indicates that the Wnt-signalling pathway is genetically buffered^^ 
and the associated genes are possibly co-evolving. 

Biological processes affected by bar-l(ga80). To investigate which 
processes were affected by the bar-l(ga80) mutation, we tested 
enrichment of mutation- affected genes in GO-, KEGG-, Anatomy- 
, Wormbook-, Gene class- and Protein domain annotations. To 
distinguish between bar-1 and developmental effects we excluded 
all the bar-l(ga80) affected genes with a developmental effect from 
the set of genes used for enrichment analysis (Supplementary 



Table 2). The results of the complete set of bar-l(ga80) affected 
genes including those with a developmental effect can also be 
found in Supplementary Table 2. 

Genes lower expressed in bar-l(ga80) are enriched with genes 
involved in cuticle constituents (p < 1*10"^), proteolysis (p < 
1*10"^°) and the proteasome core complex (p < 1*10"^°). 
Whereas the proteolysis and proteasome core complex genes over- 
lapped, they did not overlap with the cuticle constituent genes. 
Furthermore protein degradation related enrichments were reflected 
in the multiple categories tested, implying that protein degradation/ 
turn-over might be reduced. Thus, bar-l(ga80) affects protein degra- 
dation, possibly reflecting the transition of the cell from one state into 
another. 

The group of genes expressed higher in bar-l(ga80) consisted of a 
more diverse set of genes. These genes were especially related to 
transcriptional regulation, as shown by an enrichment of the GO- 
terms regulation of transcription (p < 1*10"^), sequence specific 
DNA binding (p < 1*10"^), transcription factor activity (p < 
1*10"^) and nucleus (p < 1*10"^). Some indications were found that 
the Ras-pathway was affected because transcription factors known to 
be linked to the Ras-pathway were up-regulated, like the RAS inhi- 
bitors MDL-1^' and LIN-15B^^ The activation of the Ras pathway is 
further shown by the strong up regulation of cav-1 in the bar-l(ga80) 
mutant. 

Furthermore, also neuron- related terms were represented, as 
shown by enrichments of the GO-term axon (p < 1*10"^), synapse 
(p < 1*10"^) and in the anatomy terms where the three most sig- 
nificantly enriched groups were neuronal (p < 1*10"^°). The enrich- 
ment of these neuronal genes can point in the direction of the 
aberrant neuron migration that is observed in bar-l(ga80y^, the 
mutation might affect neuropeptide signalling. 

bar-1 (gaSO) transcription patterns suggest DAF-16 activation. 

Since enrichment in transcriptional regulation was detected, we 
used the modENCODE^°'^^ set of ChlP-seq determined binding 
sites to search for enrichment of binding sites for transcription 
factors. We found that the genes higher expressed in bar-l(ga80) 
were enriched for binding-sites of transcription factors PHA-4, 
MDL-1, ELT-3 and PQM-1 (hypergeometric test, P < 1*10""). 
These transcription factors were up-regulated in the bar-l(ga80) 
mutant, except for elt-3 (Figure 4A). Together with the enrichment 
found for the binding sites, this indicates that the absence of the P- 
catenin BAR-1 results in an activation of transcription factors, 
possibly as a compensatory response. The four transcription 
factors for which enrichments have been found share binding sites 
for many of the genes. Over 50% of the up -regulated genes in the 
transcriptional network were associated with more than one of these 
four transcription factors (Figure 4B). Furthermore, PQM-1 and 
MDL-1 also bind near the transcription starting site of PHA-4 and 
ELT-3 (Supplementary figure D)^°. This indicates that it is likely that 
PQM-1 or MDL-1 is involved in the transcriptional activation 
observed in bar-l(ga80). 

Three of these transcription factors: PQM-P^ MDL-1 ^^ and ELT- 
3^^, have been associated with the insulin/IGF- 1 signalling pathway 
and longevity. However, for ELT-3 this relation is debated in more 
recent literature^^. Furthermore, PQM-1 is also identified as a pro- 
moter of growth, development and reproduction^^. PQM-1 has an 
antagonistic interaction with DAF-16, where nuclear translocation 
of PQM-1 (promoted by DAF-2) results in depletion of DAF-16 
from the nucleus (and vice-versa). Tepper et al. identified genes 
regulated by PQM-1 (referred to as class II genes), and genes regu- 
lated by DAF-16 (referred to as class I genes). It is also shown that 
some of the class I genes are also regulated by PQM-1 We tested 
expression of these genes in the bar-l(ga80) versus N2 and found 
that the class I genes were enriched for in the up-regulated genes 
(hypergeometric test, p < 1*10"^^) and slightly but significantly 
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Figure 4 | Transcription factor activity and targets in up-regulated genes. 

(A) shows the transcript abundance of the four transcription factors 
enriched for targets among the up-regulated genes in bar-l(ga80). The 
fold-change in bar-l/N2 is shown. The levels oipqm-1, pha-4, and mdl-1 
are significantly higher in bar-l(ga80) (linear model, p < 0.01), whereas 
this is not the case for elt-3 (linear model, p = 0.158). (B) A Venn-diagram 
of the up-regulated genes associated with the four enriched transcription 
factors. There is a high level of overlap between the associations as >50% of 
the targets are associated with multiple transcription factors. 

up-regulated in bar-l(ga80) (two-sided t-test, p < 1*10"^). 
Moreover, the specific PQM-l targets were enriched for in the 
down-regulated genes (hypergeometric test, p < 1*10"^) and were 
slightly down -regulated (two-sided t-test, p < 1*10"^), see also 
Supplementary figure E and Supplementary Table 3. We also ana- 
lysed DAF-16 ChlP-seq data^^ and found that the genes up-regu- 
lated in bar-l(ga80) were enriched for DAF-16 targeted genes (124 
out of 425 up-regulated genes, hypergeometric test, p < 1*10"^). 
Based on these results we hypothesize that loss of function of bar- 
l(ga80) leads to induction of a stress response reflecting DAF-16 
activation, causing delayed development of the worms. 

Conclusion 

We studied gene expression patterns during the fourth larval stage 
(L4) of the strain EW15 carrying the P-catenin-loss-of-function 
point mutation bar-l(ga80) causing a Glu to Stop codon change at 
amino acid 97 of the predicted BAR-1 protein. To untangle the 
developmental effects from the effects of the bar-1 mutation, we used 
a time-series dataset^^ We showed that bar-l(ga80) results in a 
slower development, as these worms take on average —10% more 
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^ Class I genes 
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development 



Figure 5 | A model for the bar-llof effects. This model incorporates the 
findings in this paper (grey) with what is known about the Wnt-pathway 
regarding to BAR-1 (black)^ and findings reported about DAF-16 and 
PQM-l (red)^^. It is proposed that (transcriptional) activity of the p- 
catenin bar-1 results in a feedback loop, de-activating the Wnt-pathway. 
Furthermore BAR- 1 activity is needed for a correct developmental 
program, where bar-llof shows indications of a DAF-16-mediated stress 
response. The exact level of interaction between BAR-1 and the insulin 
pathway remains to be elucidated. 



time to develop than Bristol N2. Using the transcript ome to estimate 
the age of the worms, we found that 48 hours after synchronization, 
bar-1 (ga80) worms are transcriptionally most similar to N2 worms at 
44-45 hours after synchronization^^. By measuring the time that 
egg-deposit starts in bar-l(ga80) and Bristol N2, we confirmed this 
developmental delay. To our knowledge, this has not yet been 
reported for bar-l. 

Analysis of the Wnt-pathway showed that a non-functional bar-1 
causes up-regulation of Wnt- signalling components, mom-2, cfz-2, 
mig-1, lin-18, dsh-1, mom-5 and lin-17. Together with the down- 
regulation of Wnt-inhibitor sfrp-1 this indicates hyper- activation 
of the Wnt-signalling pathway, suggesting a compensatory mech- 
anism (Figure 5). This is further shown by the modest up-regulation 
of the other P-catenins, wrm-1, hmp-2 and sys-1 which all have the 
potential to substitute for bar-1 in transcriptional activation 

Analysis of the genes affected by bar-l(ga80) showed that genes 
up-regulated in bar-l(ga80) are enriched for transcription factor- as 
well as histone- binding sites and for processes like chromosome 
rearrangement, chromatin factors and neurogenesis. The down- 
regulated genes were enriched for cuticle components and hh- 
signalling pathway genes, suggesting bar-1 directly affects these 
processes. We found that the transcriptional response induced in 
bar-l(ga80) reflects DAF-16 activation (Figure 5). This also corre- 
sponds with the developmental delay we measured. We propose that 
loss of bar-1 results in a compensatory/feedback response on the 
transcriptional level, leading to Wnt-pathway and DAF-16 activation. 
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